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Using the Keldysh nonequilibrium Green function technique, we study the current and shot noise 
spectroscopy of a single interacting quantum dot coupled to two ferromagnetic leads with different 
polarizations. The polarizations of leads can be both parallel and antiparallel alignments. General 
formulas of current and shot noise are obtained, which can be applied in both the parallel and 
antiparallel alignment cases. We show that for large polarization value, the differential conductance 
and shot noise are completely diferent for spin up and spin down configurations in the parallel 
alignment case. However, the differential conductance and shot noise have the similar properties 
for parallel alignment case in the small polarization value and for antiparallel alignment case in any 
polarization value. 
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Spin-related phenomena in semiconductor quantum 
dot have attracted great interest recently as they are the. 
crucial ingredient in the emerging field of spintronicsEJ 
and several quantum computation schemea. In addition 
to their potential industrial applications, these devices 
also provide an ideal test ground for the study of basic 
physLes including many-body effect, such as the Kondo 
cffcctB. In this Letter, we use the Keldysh nonequilibrium 
Green function technique to study the current and shot 
noise through an interacting dot coupled to two ferromag- 
netic (FM) leads as a function of the applied bias voltage 
for parallel (P) and antiparallel (AP) lead-polarization 
alignments. Using the cquation-of-motion approach, we 
obtain a general formula of current and shot noise for in- 
teracting quantum dot, which can be applied in studying 
the transport phenomena of dot coupled to FM leads with 
both the P and the AP alignments. Our results show that 
both the differential conductance and shot noise show the 
completely different behavior for spin up and spin down 
configurations in P alignment case with large polarization 
value. However, the differential conductance and shot 
noise show similar behavior for P alignment in small po- 
larization value and for AP alignment in any polarization 
value. 

The system Hamiltonian is written as 

H = H l + Hr + H d + H t ■ (1) 

The Hamiltonian for electrons in the left and right non- 
interacting metallic leads is 

H L +Hr= e kac\ a C ka , (2) 

k€L,R;cr 

where the electron creation (annihilation) operators in 
the leads are denoted by c\ a (c^). The Hamiltonian of 
the dot is 

H D = y^Jedadldg + —n d ,<xnd,a] , (3) 



where d) a (d a ) are the creation (annihilation) operators 
of dot electrons, and tda is the resonance level of the dot 
which can be tuned by magnetic field. The coupling of 
the dot to the leads is 

h t = [v4X'+^i. (4) 

where the tunneling matrix elements Vka.a transfer elec- 
trons through an insulating barrier out of the dot. 

Usinftrthe Keldysh nonequilibrium Green-.£unction for- 
malismtffl, the terminal current is given byofl 

1 = de{Tr[(f L r L - f R T R )(G r - G a )} 

+Tr[(r i -r i? ')G < ]}, (5) 

and the spectral density of shot noise in the zero- 
frequency is given byu 

S(u> - 0) = £fde{-f L (l f L )(Tr[(T L Gn 2 } 
+Tr[(r L G a ) 2 ]) + i/iTr[r L G > ] - i(l - / L )Tr[r L G<] 
+/ L Tr[r i G > r L (G r - G a )] + Tr[r i G > r L G < ] 

-(1 - / L )Tr[r L (G r - G a )r i G < ]} , (6) 

where Jl(r) & r e the Fermi distribution function of the left 
and right leads, which has different chemical potential 
upon a voltage bias \it, — fiR = eV. The coupling of the 
dot to the leads is characterized by the parameter 

T L Jf = 2^Y J PL(R)^{e)V;n >a {e)V^^{e) (7) 

a" 

with PL(R),a the spin-cr band density of states in the 

two leads. G r< " a ) and G < ^f' > are the Fourier transform 
of the dot electron retarded (advanced) Green function 

&£){t,U) - T i0(±t T i')<K(t),4'(O}). thc lesser 
Green function G< a ,(t,t') = i{dl,{t')d a {t)), and the 



2 



greater Green function G> a ,(t,t') = — i(d a {t)d [ a , it')). It 
is noted that Eqs. (||) and @ expresses the current and 
the fluctuations of current through the quantum dot, an 
interacting region, in terms of the distribution functions 
in the leads and local properties of the quantum dot, such 
as the occupation and density of states. 

In order to compute the current and the shot noise, one 
has to compute the dot electron retarded Green function 
G r and Keldysh Green function G < in the presence of 
Coulomb interaction U . Without loss of physics we are 
considering here, we assume that the tunneling matrix 
elements are spin independent, Vka-.a' = VkS aa >. Using 
the cquation-of-motion approach, we obtain the retarded 
dot Green function in the large U limit as: 



1 - {nd.a) 



J l<7 



i0+ 



(8) 



Here 



v \v k \ 2 
de r£( £ ) + r«( £ ) 



27T CO 



i0+ 



(9) 



k£L,R 



\Vk\ 2 fL/n{£ks) 



i/2r s 



de T^e)f L (e)+T R (e)f R (e) 



2tt uj - e do 



(10) 



where T^ R \e) = r^i- (e), and the occupation num- 
ber is subject to the self-consistency condition (rider) = 
-i J |f G<{lo). The finite life-time in Eq. © for finite 
bias voltage and magnetic field can be-pbtained by using 
the second-order perturbation theoryEl 

The Green function cannot be obtained by di- 

rectly using the above equation-of-motion approach with- 
out introducing additional assumptions. By instead-Ap- 
plying the operational rules as given by LangrcthLl to 
the Dyson equation for the contour-ordered Green func- 
tion, one can show the following Keldysh equation for the 
lesser and greater functions!] 



G«»(w) = C»£< (>) (w)G a M . 
From Eq. (@), the retarded self-energy is given by 
'{n ds )(u - e da ) + (YjQ a + E r lcr ) 



r(") 



We then arrive at 



1 - {nds) 



(11) 



(12) 



= i(l-(n d s))- i m^)+T^)) 



-2ImE[ CT H} 



(13) 



To determine and , we follow the ansatz pro- 
posed in RefO to assume further that these self-energies 
have the form 

S< = i[T L a {io)f L (u) + r?(w)Mw)]R„(u,) , (14a) 



= -i[Tj(a>)(l - h{u)) + r?(w)(l - /fl(w))]i^(a;) , 

(14b) 

where R(u>) can be regarded as a renormalization factor 
due to strong electron-electron Coulomb repulsion in the 
dot. A little algebra yields 

= i _ L-\ ^ - r7 ImEL (a;)} (15) 

where T ff = + r£. 

As an extension, we can assume in general the "scat- 
tering in" and "scattering out" self-energies to have the 
form 

£< = i[T L (u>)f L (u}) + T R {w) f R (u>)]TL{w) , (16a) 



£> = -i[T L (u;)(l - f L (uj)) + T R (u;)(l - f R (u))]R{ui) , 



Eqs. (16a) and (16b) leads to 



G r - G a = G> G' 



(16b) 



(17) 



Substitution of Eq. (|17|) and Eq. (|llj) with Eq. (|16|) into 
Eqs. (0) and m) gives rise to 



-iG r (T L +r H )RG° 



I 



± J ATr{r L [G<+/ L (G' -G a )]} 

± J deTr{T L [iG r (f L T L + f R T R )KG a 
-if L G r {T L +r K )RG Q ]} 
^ | - /«]Tr[T] , 



(18) 



and 



S = £f &{[/L(l-/L)Tr{T[l + 2(r^R)- 1 (R-l)]} 
+Ir(1 - f R )Tr[T] + (f L - f R ) 2 Tr[(l - T)T]} (19) 

where the transmission coefficient matrix is defined as 
T = G c T L G r r if R. Wc remark that Eqs. @ and 
( |l9| ) are the general formulas for current and shot noise 
in the interacting dots with large U, which can be ap- 
plied in studying the spin-transport of both the P and 
AP configurations. It is noted that the formula of pro- 
portional coupling (i. e., = XT R ) does not work in 
the AP alignment case at finite bias voltage and mag- 
netic field. However, one can easily to compute the cur- 
rent and shot noise for AP alignment configuration by 
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FIG. 1: The zero-temperature differential conductance of spin 
up (red-solid line) and spin down (blue-dashed line) for FM 
leads with P alignment as a function of bias voltage without 
(a) and with (b-f) Zeeman splitting Ae^ = e^t — e <ii- Here 
ro = 1, the dot level in the absence of magnetic field is ej-f = 
edi = —3.0, and the band width is 100. Insets in (b,d,f) 
indicate the DOS of spin down configuration as a function of 
bias voltage near eV — 0. 



using our formalisms Eqs. @ and @. It is easy 
to find that for the noninteracting case, the above two 
expressions are just the Landaucr-Buttikcr formalisms 
developed based on the scattering matrix theory. The 
connection between the two formalisms was fkst estab- 
lished by Mcir and Wingreen for the currentB, and by 
Zhu and Balatsky for the shot noisecl. 

In the wideband limit, we assume that Tv {R} to be 
energy independent, but polarization dependent with 
(1 + P)r , rf = Tf = (l - P)T for P align- 



1 1 



1 T 



mcnt, while T 



(i + p)r , rf = rf = (i-P)r 



for AP alignment, where Tq describes the tunneling cou- 
pling between the dot and the nonmagnetic leads, and 

< p < 1 characterizes the polarization of leads. In Fig. 

1 we plot the zero-temperature differential conductance 
of spin up (solid line) and spin down (dashed line) for FM 
leads with P alignment as a function of bias voltage. One 
interesting observation is that the behavior of differential 
conductance is quite different for small and large P val- 
ues. For small P = 0.2, Fig. 1(a) shows that the zero- bias 
peak of differential conductance for nonmagnetic leadstHl 
splittes into two peaks even in the absence of magnetic 
field, and this splitting can be tuned away by applying an 
appropriate magnetic field [Fig. 1(c)] and can be recov- 
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FIG. 2: The differential conductance for FM leads with AP 
alignment. 



ered by increasing the magnetic field furthcr_JFig. 1(c)]. 
This result agrees well with that in Ref. E3, which is 
due to the splitting of the dot levels renormalized by the 
spin-dependent interacting self-energy Zdcr satisfying the 
self-consistent equation Ida = e^ + ReE^ (eda,£da)- For 
large P = 0.8 [Figs. 2 (b,d,f)], the behavior of spin up is 
similar to that of small P = 0.2, which shows one peak 
at approriate value of magnetic field and two peaks when 
the field is away from the appropriate value. However for 
the spin down configuration, the differential conductance 
shows a small "flat" near eV = when the magnetic 
field is smaller than the appropriate value, while it shows 
a "dip" at eV = when the field is at the appropriate 
value, and finally it shows a large maximum when the 
field increases further. This interesting behavior can be 
understood with the help of DOS for spin down near the 
zero bias voltage (shown as insects of Figs. l(b,d,f)). 
The DOS for spin down configuration shows a dip at 
lo = €di — ?dT because Re E^, (u>) grows logarithmically 
at this value, and the dip position can be tuned by mag- 
netic field around the appropriate value corresponding to 
the single peak for spin up configuration. In Fig. 2 we 
plot the differential conductance for FM leads with AP 
alignment. In this case, the behavior of differential con- 
ductance looks similar for spin up and spin down config- 
urations, which has one peak in the absence of magnetic 
field [Fig. 2(a,b)] and two shifted peaks in the presence 
of magnetic field [Fig. 2(c,d)]. The small difference be- 
tween the up and down spin configurations is due to the 
small (compared to the P alignment case) splitting of 
renormalized dot levels for different spin configurations 
at finite voltage (as shown in the inset of Fig. 2(b)). 

In Figs. 3 and 4, we plot the corresponding differential 
shot noise for FM leads with P and AP alignments as a 
function of the bias voltage. For P alignment case, the 
differential shot noise Fiq. 3 also shows the different be- 
havior for spin up and down configurations in large P val- 
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FIG. 3: The differential shot noise for FM leads with P align- 
ment. 
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FIG. 4: The differential shot noise for FM leads with AP 
alignment. 



ues, while similar behavior for small P value. However, 
for AP alignment case, the differential shot noise Fig. 4 
shows the similar behavior for different spin configura- 
tions. Since the transmission probability corresponding 
to the conductance peaks (as shown in Figs. 1 and 2) is 
small, the differential noise is approximately proportional 
to the conductance and only shown one single peak rather 
than two peaks. 



To summarize, using the Keldysh nonequilibrium 
Green function technique, we have studied the current 
and shot noise spectroscopy of a single dot with Coulomb 
interaction coupled to FM leads with P and AP polar- 
ization alignments. We have shown that the lead align- 
ments affect both the current and current fluctuations. 
For large polarization value, the spin up and spin down 
configurations have the different behavior in the differen- 
tial conductance and shot noise as a function of bias volt- 
age in P alignment case. While the differential conduc- 
tance and shot noise show similar behavior for different 
spin configurations in P alignment with small polariza- 
tion value and in AP alignment case with any polariza- 
tion value. The derived current and shot noise formulism 
can be applied to more complicated system. Work along 
this line is still in progress. 
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